Mesoscopic charging effects in the a-model for granular metals 
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We derive the a-model for granular metals. It is valid for any relation between the temperature, 
T, and the grain mean level spacing. The continuum limit of the cr-model describes non-perturbative 
charging effects in homogeneous disordered metals. The cr-model action contains a novel term which 
is crucial for studying Coulomb blockade effects. Topological properties of the cr-model encoding 
charge discreteness are studied. For T below the escape rate from the grain, F, Coulomb blockade 
effects are described by solitons of the Q-matrix which are delocalized in real space. At T > F the 
solitons transform into phase instantons localized at a single tunneling contact. 

PACS numbers: 73.23Hk, 73.40Gk, 73.22Lp 
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Physics of electron-electron interactions in mesoscopic 
systems has been the subject of intensive experimen- 
tal and theoretical studies 0. Recently, much at- 
tention was attracted to the study of granular met- 
als H H 0, H SH . They can be viewed as extended 
networks of metallic grains connected by tunneling con- 
tacts with dimensionless conductance qt- The Coulomb 
interaction in such systems may be described by capaci- 
tive coupling between the grains, 



He = 



1,1' 



(1) 



Here e is the electron charge, 1 and 1' label the grains, 
Cw is the capacitance matrix, iVi is the electron number 
operator in the 1-th grain, and qi are linear in external 
electrostatic potentials, such as gate voltages. 

At temperatures below the charging energy, Ec ^ 
e^/C, the fluctuations of the number of electrons in the 
grain are suppressed and electron transport is inhibited 
due to Coulomb blockade effects. The latter are typ- 
ically treated d, IM using the dissipative action of 
Ref. llol| in terms of the phase 0(r) conjugate to the grain 
charge. Coulomb blockade effects are non-perturbative in 
the tunneling conductance gx and are described by the 
phase instantons The approach of Ref. neglects 
elastic propagation of electrons across the grain (elas- 
tic cotunnelin g) a nd is valid for zero mean level spacing, 
i5 — > 0, 113L I14II . At temperatures smaller than the es- 
cape rate from the grain, F ~ gr^, elastic propagation of 
electrons becomes imp ortant. In this regime the dissipa- 
tive action of Ref. should be replaced by the cr-model 
in which electron propagation across grains is described 
by the Q-matrix. 

For gx ^ 1 the electron transport at T <C F is diffu- 
sive . Granular metals in this regime may be described 
by the continuum limit of the tr-model. It is expected 
that the continuum limit of the a-model is independent 
of the granular nature of the metal. Therefore, study- 
ing charging effects in this regime will shed light on the 
role of charge discreteness effects in homogeneous disor- 



dered metals which are expected to become important 
close to the metal-insulator transition. In the metallic 
regime Coulomb blockade effects are non-perturbative in 
l/gr and in the cr-model treatment should be described 
by Q-matrix solitons which replace the phase instantons 
at T < F. An important step in this direction was made 
in Ref. . Since the consideration of Ref. 0| assumes 
zero mean level spacing in the grain and spatial sepa- 
ration between the Coulomb interaction and Q-matrix 
degrees of freedom it is not clear how to generalize it to 
the continuum limit. 

Here we derive the a-model for granular metals, 
Eq. l(TT|l . It enables one to study non-perturbative charg- 
ing effects at T < F. In contrast to Refs. 0, Q its ap- 
plicability is not limited by the condition T ^ 6. The 
continuum limit of the cr-model, Eqs. (|15|l and H16() . de- 
scribes charging effects in homogeneous disordered met- 
als. Equations IjlSI) and represent our main 
results. 

We consider a d-dimensional array of metallic grains 
coupled by tunneling contacts with dimensionless con- 
ductance gx — , with R being the contact resistance. 
We assume that the Thouless energy of the grain, Et^ 
exceeds both Ec and the elastic escape rate, F. We 
make no assumptions about the relation between T and 
6, where S is the mean level spacing of the grain per 
spin projection. For simplicity of presentation we re- 
strict ourselves to the case of normal metals with bro- 
ken time reversal symmetry. The generalization to other 
symmetry classes and superconducting case will be con- 
sidered elsewhere. We consider electrons moving in the 
presence of a random potential, Vi(r), with a variance 
(Fi(r)Vi' (r')) = 2^S{r — r')6\i', where v is the density 
of states and r is the mean free time. Electron tunnel- 
ing between adjacent grains is described by the tunneling 
Hamiltonian 



l\',pk 



(2) 



where cj^ is the electron creation operator. Finally, the 
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Coulomb interaction is described by Eq. 

Thermodynamic and linear response transport proper- 
ties of a granular metal can be obtained from the aver- 
age replicated partition function, {Z"-) = ^exp a-^^^, 

where a is the number of replicas and the averaging is 
performed over the random impurity potential. We write 
(Z") as an imaginary time functional integral over the 
fermions. Averaging over disorder induces a quartic in- 
teraction of fermions which is decoupled using a Hcrmi- 
tian matrix field Q [ifil fl^ . The Coulomb interaction 
is decoupled by the scalar field, in each replica. 

Integrating out the fermions we obtain. 



Tc = 



/3 



ill' □ 



^V^ji' -I- ibwqx 



-Trln {ii^iY-i-H^^Y 



(3a) 
(3b) 

(3c) 



Here e = id^, and -ffo — ~ with /i being the 
chemical potential, is the single particle Hamiltonian in 
the absence of tunneling and disorder. The operator t 
describes tunneling between grains, and we omitted the 
argument of the scalar potential Vji(t). The trace in 
Eq. (j5c|l is taken over the spin, replica, Matsubara and 
grain indices. 

It is convenient to introduce the phases, Xii(''')i by sin- 
gling out the static part, V^^, of the auxiliary field. 



^.iM=^°+X,i(t), 



(4) 



where Xii(O) = Xji(/5)- The static part, V'^^, plays the role 
of fluctuating chemical potential and is responsible for 
charge neutrality. The time-dependent part of Vj\ (t) can 
be removed from Eq. (|3cl) by the gauge transformation. 



g(r,r') 



exp[zx(T)]Q(T, T')exp[ 
iexp[-ixjii' (t)]. 



'*x(t')], 



(5a) 
(5b) 



where xu' denotes a diagonal matrix of phase differences, 

Xii' {t) = Sij [Xji (t) - Xji' (t)] ■ I— 1 I— 1 

Following the standard procedure |l6l [iTj we look for 
diagonal coordinate-independent saddle points of the ac- 
tion in Eq. H3a(l . The tunneling Hamiltonian in Eq. 13c() 
leads to a finite escape rate from the grain, T ^ t'^ /S. For 
r <IC l/r this term may be ignored while looking for the 
saddle point. 



Qj\ = 



li + iVj'i - Ho 



iQj\ 
2t 



(6) 



We denote by A the solution of Eq. © which for En <^ fJ. 
has the form, (5ji(e„) = A(e„) — sign(e„). It has the 



meaning of the average equal point Green's function of 
an isolated grain with the chemical potential fj, + iV^^ . 

At high energies, |e„| ^ 1/r the fluctuations of the Q- 
matrix about the saddle point are strongly suppressed, 
and it can be replaced by Q = A. Using this we separate 
the high energy part of the action, Eq. ll^ . and rewrite 
the effective low-energy action in the form, 



MV] - ^^Tr[e**"'A,e-^*"'Ap 
11' 



-Tr' In 



■ + - i 



ie + tV^ -t- Ho + i^ 




(7) 



To arrive at Eq. |7J) we expanded the right hand side of 
Eq. to second order in t. Here Tr' denotes the trace 
restricted to the low-energy sector, |e„| < l/r, and .7^o[^] 
denotes the action Td, Eq. iPcJ, in the absence of tun- 
neling, t — 0. We may expand ^o[^] to second order in 
Taking into account that the electron charge den- 
sity is cancelled by that the ionic background we obtain. 

The action, Eq. |(7J), is minimized at the saddle point 
given by the solution of Eq. 10. The most general di- 
agonal solution is given by Qo{en) = il- The massive 
fluctuations about the saddle point can be integrated out 
in the Gaussian approximation. The remaining integral 
over the massless modes constitutes the cr-model. For 
energies e ^ Et the massless modes are parameterized 
by Q-matrices which are coordinate independent within 
each grain (zero-mode approximation). They can be ex- 
pressed as 



Qi - uIQoUu 



(8) 



where Qq is a diagonal matrix with diagonal elements 
equal to ±1, and Ui are unitary matrices which differ 
from unity only for |e|, |e'| < Et- Thus the manifold of 
massless modes breaks up into disjoint manifolds charac- 
terized by the trace of the Q-matrix, 



TrQi = 214^1. 



(9) 



Restricting ourselves to the massless mode manifold, 
Eq. ijSJ), we expand the last term in Eq. (TJ to linear 
order in e and and to second order in t, 

Ti = To[0] - ^ETr[e*"Qie"*"Qi'] 



11' 



7E(^ + Tr'[(e + V^^)/^.(Qi-A)]), (10) 



0\2 



ttT 



where Pj is the projection operator onto the j-th replica. 
Substituting Eq. H10|) into Eq. integrating over V^^ 
and taking into account that Tr'[PjA] = we obtain the 
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cr-model for a granular metal, 




Here the trace is taken over the spin, replica, and Mat- 
subara indices. The summation in Eq. (|lla|l goes over 
the integers, Wi in each grain, see Eq. ©. In Eq. ifTT)) 
introduced the notation, 

N,i = -i_Tr[Q,P,] - qi. (12) 

Equation (|12|) has a very natural interpretation: the 
Q-matrix parametrizes fluctuations of electron Green's 
function, and Nj\ represents the deviation of the charge 
of the 1-th grain in replica j from qi due to Q-matrix 
fluctuations. 

We emphasize that we did not have to assume T » 
S in order to integrate out the static component of the 
auxiliary field VR for a fixed value of the Q-matrix. The 
cr-model, Eq. (|ll|l . is valid for any relation between T and 
6. Its action differs from those in Refs. 0,01 by the second 
term in Eq. Ijllb|l . This term describes charge neutrality. 
It is instructive to show how this term manifests itself in 
the description of charging effects. 

First let us consider isolated grains, gx = 0, and set 
qj\ 0. In this case the integral over the phases x decou- 
ples from the Q-matrix integral. The action in Eq. Ill bp 
is minimized when the Q-matrix is diagonal in replica, 
Matsubara, and spin space, Qi = A (5,^ (5„m(5^i,sign(e„), 
where e„ — 7rT(2?i + 1) are fermionic Matsubara fre- 
quencies. Using the standard parametrization for the 
Q-matrix 17], Q = e™Ae~™, we observe that the ex- 
pansion the second term in Eq. (|llb|l starts with the 
fourth order in u. At T 3> (S the fluctuations of u are 
small, (u^) ^ 6/T, and the second term in Eq. Ijllb|) 
can be treated perturbatively. As a result we obtain the 
following correction to the average free energy, 

^^-i(¥)^i'-^i:((^i*^'i)^)„^ (13) 

Here (. . .)o denotes the averaging with respect to the 
zero-dimensional action given the first term in Eq. (|llb|l . 
Equation H13|l has a clear meaning. Coulomb interac- 
tion enforces charge neutrality. Thus the grain charge 
is independent of the disorder realization is determined 



by the external electrostatic potential. Therefore ensem- 
ble averaged thermodynamic properties can be replaced 
by canonical averages for non-interacting electrons. It is 
well known that for the canonical ensemble mesoscopic 
fluctuations result in the followin g co rrection to the av- 
erage free energy of the grain [13, Hfl , 

AF^{5/4)({ANr)„ (14) 

where AiV is the fluctuation of the number of parti- 
cles in the grain at a constant chemical potential /i. 
The correction ()14|l is crucial for some quantities, such 
as persistent currents 0, 0, I20I I2H The quantity 
—i^TrlQiPj] gives the fluctuation of the number of 
particles in the j-th replica in grain 1, see Eq. (|12ll . 
Therefore Eq. H13|l is equivalent to Eq. H14|l . Perturba- 
tive evaluation of Eg. ifT!^ is straightforward and gives 
the result of Ref. \M for a zero-dimensional system, 
A^^ = ~d/i27T^)JJdede'ie-er'fie)fie'), with /(e) 
being the Fermi function. 

Apart from the perturbative correction, Eq. (|13|) . the 
presence of the second term in Eq. Ijll b|) leads to topolog- 
ical consequences which are important for the description 
of Coulomb blockade effects. In the absence of tunneling 
the action in Eq. (jllb|l has a discrete set of degener- 
ate minima which are achieved at diagonal Q-matrices, 
(3i(e) = SijSnmSp,i,sign{en — 2'nTwj\), where Wji are inte- 
gers analogous to the winding numbers in the AES ap- 
proach. It is easy to check the degeneracy of these min- 
ima as the dependencies on Wji in the first and second 
terms in Eq. Ijllb|) exactly cancel each other. This is a 
straightforward consequence of gauge invariance since the 
different minima are related by singular gauge transfor- 
mations, (3(T,r') ^ e*2'^^"'^0(T,r')e-'27rT™r'^ rpj^g g^j^ 

2 Wj^i = W\, where the factor of two arises from spin 
degeneracy, defines the trace of the Q-matrix via Eq. 
Different minima with the same W\ belong to the same 
manifold of Q-matrices and can be connected by contin- 
uous rotations in the replica and Matsubara space. 

At T ^ ^ the Q-matrix fluctuations around the diago- 
nal minima are suppressed. One can integrate over them 
in the Gaussian approximation, and replace the integral 
over Q in Eq. Hlla|) by the sum over the winding num- 
bers Wjx. This summation is equivalent to the summation 
over all charge states of the grain . li gr is finite but 
the temperature is much higher than the elastic escape 
rate from the grain, T ^ F ~ 5t<5, the tunneling term in 
Eq. fTTH|) may be considered as a perturbation. It merely 
breaks the degeneracy of the minima. In this limit one 
recovers the action of Ref. |0| from Eq. (|ll|l . Coulomb 
blockade effects in this regime are described by discon- 
tinuous changes of the winding numbers Wj^x which occur 
at the tunneling contacts between grains and correspond 
to instantons of Ref. [Il| . 

At T < F the tunneling term in Eq. I|llb|l may not 
be treated perturbatively. Since this term is sensitive 
only to the spatial variations of the Q-matrix and of the 
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phases x it preserves the degeneracy of spatially uni- 
form minima of the action in Eq. Ijllb|) : x = and 
Q = (5y 5„mSign(e„ — 2TrTwj\). The Q-matrix may fall 
into different minima in different spatial domains. If the 
sum 2 Wji — Wi in the two domains is the same 
the two minima belong to the same Q-matrix mani- 
fold. The Q-matrix configuration which minimizes the 
action, Eq. (|llb(l . at the boundary between two such do- 
mains corresponds to a gradual soliton-like rotation of 
the Q-matrix in replica and frequency space. The pre- 
cise form of such solitons depends on the dimensionality 
of the system. Their spatial extent can be estimated 
from Eq. (|irbl) to be the temperature diffusion length, 
yryr, where a is the grain size. At T F the 
solitons are delocalized over many grains and can be de- 
scribed by the continuum limit of the a-model, Eqs. I|15|) 
and (1161). Such solitons are analogous to those found in 
Ref. Ell. 

The cr-model, Eq. has a well defined continuum 
limit. It is obtained by assuming slow spatial variations 
of the Q-matrix and of x = ^ijXji expanding the last term 
in Eq. (jllb|l to second order in gradients, and replacing 
summation over 1 by integration over r = al. The tun- 
neling term in Eq. I|llb() in the continuum limit describes 
diffusion with the diffusion constant D = ^qtSo?', where 
the coefficient 7 depends on the spatial arrangement of 
the grains, 




where dQr = VQr — *[Vxr, Qr], and denotes the local 
shift of chemical potential induced by the fluctuation of 
the particle number, Njr^ 

S T S 

= = -i^TviQrP,] - -Qr. (17) 

As TiQ is an integer, see Eq. (jSJ , it must remain constant 
in space in the continuum limit. Different W's reflect 
charge quantization in the whole system. Such effects 
may be disregarded in the thermodynamic limit, and the 
Q-integral should be taken over the traceless Q-matrices. 

In conclusion, we derived the cr-model for granular 
metals, Eq. Hll|l . It describes charging effects in gran- 
ular metals in the low temperature regime which in- 
cludes T < S. The (T-model action Hll(l contains a novel 
term whose presence is crucial for describing charging 
effects. For an isolated grain perturbative treatment of 
this term reproduces the canonical ensemble result (|14|l . 
In the limit i5 — > the a-model reproduces the action 
of Ref. i'^is limit the Coulomb blockade effects 

are described by phase instantons of Ref. ^3- -^t tem- 
peratures below the escape rate from a grain, F ~ grS, 



Coulomb blockade effects within the cr-model framework 
are described b y so litons of the Q-matrix similar to those 
studied in Ref. . The spatial extent of such solitons is 
given by the temperature diffusion length, Lt = a^jT /T. 
Inclusion of such solitons may be important for studying 
charging effects in granular systems at T < F. 
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